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A  LIFTING-SURFACE  PROGRAM  OF  CONTRAROTATING  PROPELLERS 

Beniamin  Y.-H.  Chen  and  Arthur  M.  Reed 
David  Taylor  Research  Center 
Bethesda,  Maryland  20084-5000 


ABSTRACT 

A  new  lifting-surface  computer  program  for  a  set 
of  CR  propellers  has  been  developed  based  on  a  modi¬ 
fied  version  of  the  MIT  lifting-surface  design  program 
with  hub  effects.  This  program  automatically  computes 
the  velocities  Induced  by  one  propeller  on  the  other. 
In  addition,  the  hub  portion  of  the  program  is 
modified  to  account  for  the  velocities  induced  by  one 
propeller  on  Che  hub  of  the  opposite  propeller.  Data 
from  LDV  measurements  of  induced  velocities  have  been 
used  to  adjust  the  shape  and  distribution  of  the  wakes 
shed  from  the  two  propellers.  A  comparison  between 
the  conventional  and  the  new  methods  for  the  design 
of  a  sec  contrarotating  propellers  for  a  surface  ship 
is  also  given. 

I.  INTRODUCTION 

Over  the  past  several  years,  the  application  of 
contrarotating  (CR)  propellers  on  surface  ships  has 
been  reemphasized  because  of  interest  in  utilizing 
lighter  and  more  efficient  electrical  propulsion  sys¬ 
tems.  This  interest  has  been  heightened  by  the 
potential  of  CR  propellers  to  achieve  Increased 
propulsive  efficiency  and  Increased  cavitation  incep¬ 
tion  speeds.  A  CR  propeller  set  is  shown  in  Fig.  1. 

Four  decades  ago  most  research  on  CR  propellers 
was  experimental,  and  design  techniques  were 
empirical.  Since  that  time,  analysis  techniques  have 
developed  so  that  today  there  are  a  number  of  sophis¬ 
ticated  lifting-line  and  lifting-surface  techniques 
available  for  use  in  design.  One  of  the  earliest 
design  techniques  for  CR  propellers  was  developed  by 
herbs  (2).  herb's  method  was  further  refined  into  a 
usable  design  technique  by  Morgan  (3).  Both  herbs 
and  Morgan  relied  on  classical  circulation 
theory  (4, 5)  for  single  propellers  as  the  basis  of 
their  methods.  To  determine  the  forces  and  induced 
flow  field,  they  applied  lifting-line  theory  for 
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moderately  loaded,  wake  adapted  single  propellers  to 
the  forward  and  aft  propellers  individually  so  that 
the  force  and  induced  flow  field  could  be  determined. 
The  variations  of  inflow  velocity  (time-average)  in 
the  axial  and  tangential  directions  and  the  streamline 
contraction  in  the  slipstream  of  the  forward  propeller 
were  treated  approximately  to  account  for  the  mutual 
Interactions  between  the  two  propellers. 


FIC.  1.  A  CR  PROPELLER  SET 


Following  Morgan's  work,  two  computer  programs 
were  developed.  Caster  and  LaFone  ^6)  developed  a  CR 
lifting-line  program  based  on  the  lifting-line  program 
for  single  screw  propellers.  This  program  incorpor¬ 
ated  field  point  velocity  computations  to  determine 
the  Interaction  velocities  between  the  two  proximate 
propellers  in  an  iterative  procedure.  Nelson  ( 7_) 
developed  lifting-line  as  well  as  lifting-surface 
programs.  He  showed  a  successful  design  for  CR  pro¬ 
pellers  on  torpedoes  using  the  circulation  distribution 
with  finite  values  at  the  root.  Both  Caster's  and 
Nelson's  lifting  life  programs  use  the  same  approach, 
but  differ  somewhat  in  the  options  that  are  available 
to  the  user  and  the  techniques  by  which  the  details 
are  carried  out. 

Cox  (8)  chose  the  Cascer  and  LaFone  computer  code 
as  the  basis  of  his  method,  but  with  several 
revisions  as  follows:  1  -  Cox  used  the  induced  axial 
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velocities  In  the  propeller  planes  to  determine  the 
mass  flow  rate  for  each  propeller,  and  allowed  that 
the  diameter  of  the  downstream  propeller  to  be  speci¬ 
fied  directly,  regardless  of  the  mass  flow  rate. 
Caster’s  program  used  an  approximation  due  to  Lerbs  to 
determine  the  appropriate  diameter  for  the  aft 
propeller;  2  -  Both  programs  use  FPV-7  to  calculate 
the  velocities  induced  by  one  propeller  on  the  other. 
However,  Cox  created  these  velocities  at  radii  smaller 
than  the  diameter  of  the  forward  propeller  hub.  Thus 
his  program  could  readily  deal  with  tapered  hubs, 
which  Caster's  program  could  not. 

Another,  more  recent,  lifting-line  computer 
program  for  the  design  of  CR  propellers  was  developed 
by  Reed  (9).  It  was  also  baaed  on  the  Caster  and 
Lafone  method,  and  it  incorporated  many  of  the 
improvements  due  to  Cox.  In  addition,  the  program 
Incorporated  Improved  methods  for  computing  the 
velocities  induced  by  one  propeller  on  the  other. 

These  new  methods  have  decreased  the  time  required  for 
the  induced  velocity  calculations  by  more  than 
90  percent,  substantially  reducing  the  time  required 
for  the  total  lifting-line  design. 

In  the  preliminary  design  stage,  lifting-line 
theory  is  a  very  useful  tool.  The  dependence  of 
forces  and  efficiency  on  such  parameters  as  angular 
velocity  diameter,  blade  outline,  blade  numbers,  and 
circulation  distribution  can  be  determined  economi¬ 
cally  through  parametric  calculations.  Since  CR 
propeller  design  involves  a  large  number  of  design 
parameters  and  the  calculation  of  mutual  interactions 
between  two  propellers,  an  economical  lifting-line 
tool  is  particularly  important  for  CR  propeller  design. 
Lifting-line  theory  alone,  however,  can  neither  deter¬ 
mine  the  final  meanline  distribution  and  radial  pitch 
variation  accurately  nor  can  it  predict  interactions 
between  two  adjacent  propellers  satisfactorily.  For 
the  final  design  stage,  lifting-surface  theory  must 
be  employed  in  order  to  incorporate  three-dimensional 
flow  field  effects.  The  traditional  design  procedure, 
described  in  detail  in  section  II,  employs  a  single 
propeller  lifting-surface  program  to  determine  blade 
final  pitch  and  camber  for  each  blade  row  separately. 

In  this  traditional  procedure,  the  interaction 
velocities  between  the  propellers  are  those  from  the 
lifting-line  model,  not  the  Lifting-surface  model. 

Also,  this  procedure  will  not  allow  zero  net 
circulation  at  the  hub  for  the  two  propellers. 

More  recently,  a  customized  lifting-surface 
computer  program  for  the  detailed  design  of  the 
individual  propellers  of  a  CR  set  has  been  developed 
to  strengthen  the  traditional  design  procedure.  This 
design  tool  automatically  computes  the  velocities 
induced  by  one  propeller  on  the  other.  Additionally, 
it  accounts  for  the  hub  of  each  propeller  in  the 
design  of  the  set  (a  factor  which  is  of  greater 
Importance  for  CR  propellers  chan  for  single 
rotation  propellers),  and  allows  finite  circulation 
at  the  hub  of  each  propeller  (the  net  circulation  for 
the  two  propellers  at  the  hub  should  be  zero).  The 
shape  and  distribution  of  the  wakes  shed  from  Che  two 
propellers,  which  is  particularly  important  for  Che 
aft  propeller,  is  adjusted  using  data  from  Laser 
Doppler  Velocimetry  (LDV)  measurements  of  induced 
velocities  around  a  pair  of  CR  propellers. 

A  sample  calculation  shows  the  comparisons  between 
the  conventional  and  the  new  method  on  a  CR  propeller 
design  for  a  surface  ship. 

II.  METHODOLCY 

As  stated  before,  one  of  the  major  casks  in  CR 
propeller  design  is  to  obtain  the  induced  velocities. 


which  include  the  9elf-induced  velocity  and  the 
induced  velocity  from  the  other  blade  row.  Fig.  2 
shows  the  velocity  diagram  for  a  CR  propeller  set. 
The  details  of  the  traditional  design  procedure  and 
the  new  design  procedure  will  be  described  in  this 
section. 


Forward  propeller  (right-handed) 


Aft  propeller  ( lef t-harded ) 


FIG.  2  VELOCITY  DIAGRAM  OF  A  CR  PROPELLER  SET 


1.  Traditional  Design  Procedure 

A  flow  chart  for  the  traditional  design  procedure 
for  a  CR  propeller  set  is  shown  in  Figs.  3  and  U.  In 
Fig.  3,  the  lifting-line  program  is  used  to  design  the 
forward  propeller  using  the  inflow  from  wake  measure¬ 
ments.  Without  considering  the  effect  of  the  aft 
propeller,  the  results  of  the  lifting-line  program  are 
the  characteristics  of  the  forward  propeller.  Mean¬ 
while,  the  field  point  velocities  Induced  by  the 
forward  propeller  on  the  aft  propeller  are  computed  by 
the  internal  field  point  velocity  program.  The 
diameter  of  the  aft  propeller  is  determined  based  on 
mass  flow  conservation,  and  the  characteristics  of  the 
aft  propeller  are  obtained  through  lifting-line  calcu¬ 
lations  using  the  combination  of  the  induced  velocities 
from  the  forward  propeller  and  the  wake  measured  at 
the  aft  propeller  position.  Once  the  aft  propeller 
lifting-line  calculation  is  complete,  the  field  point 
velocities  induced  by  the  aft  propeller  on  the  forward 
propeller  are  computed  and  Incorporated  into  the  wake 
for  the  forward  propeller.  The  lifting-line  design  is 
iterated  through  the  above  procedure  a  second  time  to 
account  for  the  interaction  effects  between  both 
prope 1 lers . 
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FIG.  3  FLOW  CHARI  uF  CZ  LIFTING-LINE  PROCRAM 


•FIG.  4  FLOW  CHART  OF  CONVENTIONAL  CR 
LIFTING-SURFACE  PROGRAM 


As  shown  in  Fig.  4,  the  resales  of  Che  lifting¬ 
line  calculations  provide  information  (principally 
circulation,  hydrodynamic  pitch  distribution,  and 
induced  velocities  between  blade  rows)  to  the  single 
propeller  lifting-surface  program  which  determines 
the  f<nal  blade  pitch  and  camber.  Two  questions  are 
raiseu  by  using  the  above  design  procedure.  First, 
the  onset  flows  of  the  lifting-surface  procedures  are 
covered  over  from  the  lifting-line  procedure.  In 
other  words,  the  field  point  velocities  induced  by  one 
propeller  on  the  other  are  based  on  the  preliminary 
blade  characteristics  through  lifting-line  design. 
Gecond,  che  slipstream  contraction  of  a  CR  set  should 
be  different  from  that  of  a  single  propeller.  To 
improve  upon  the  above  procedure,  a  new  design 
methodology,  described  In  this  report,  has  been  devel¬ 
oped  . 


2.  New  Design  Procedure 

The  CR  lifting-surface  program  Is  based  on  the 
MIT  lifting-surface  computer  program,  PB^-11  (10) 
which  includes  huo  image  effects.  PBD-11  deteTmines 
che  propeller  blade  shape  for  a  prescribed  circulation 
distribution  and  a  given  hub  geometry.  The  vortex 
lattice  approach  is  used  to  represent  the  blades  and 
their  wakes.  The  hub  is  represented  by  a  distribution 
of  vortices  which  ends  at  the  hub  apex.  By  taking 
account  of  the  hub  effect,  PBD-11  allows  non-zero 
circulation  at  the  hub  of  a  propeller,  a  factor  which 
seems  to  be  more  significant  for  CR  propellers  than 
for  single  rotation  propellers.  Additionally,  its 
wake  model  is  more  realistic. 

A  vortex/source  lattice  method  has  been  employed 
in  the  present  field  point  velocity  computation 
scheme  (Kerwin  (1_1)).  The  discretized  version  of  a 
propeller  Is  comprised  of  a  blade,  Its  wake  and  the 
hub.  Following  Greeley  and  Kerwin  (^2_),  the  wake  is 
composed  of  a  transition  wake  and  an  ultimate  wake. 

The  vortex  sheet  tends  to  contract  and  roll  up  in  the 
transition  wake.  A  single  helical  tip  vortex  and  a 
hub  vortex  forms  the  ultimate  wake. 

The  details  of  the  vortex/source  lattice  method 
are  described  as  follows.  The  line  source  strength, 
qn  m  due  to  thin  wing  theory  and  the  spanwise  vortex 
strengths,  ^m(s),due  to  the  normal  boundary  condi¬ 
tion  constitute  the  primary  singularities.  The  first 
index,  n,  stands  for  chordwlse  position  and  the  second 
index,  m,  spanwise  position.  Based  on  the  conservation 
of  vorticity,  the  secondary  singularities  include  the 
chordwise  vortices,  rnm(c),  the  transition  wake 
trailing  vortices,  rnm(tw),  the  ultimate  tip  vortex, 
rn(t^  and  the  hub  vortex,  ^(h),  where  j  is  the 
jth  panel.  A  special  treatment,  shown  in  Greeley  and 
Kerwin  (12)  and  in  Kerwin  and  Lee  (1^),  is  required 
for  the  chordwise  vortices  originating  from  the  outer 
end  of  the  tip  panel. 

Induced  velocities  due  to  the  blade,  its  vortex 
wake  and  the  hub  can  be  found  by  summing  the  product 
of  the  singularity  strengths  with  the  corresponding 
velocity  influence  functions,  H^m  ‘'and  Velocity 

influence  functions  are  defined  as  vector  velocities 
which  a  line  vortex  and  source  of  unit  strength  induce 
at  the  field  point.  These  Influence  function  veloci¬ 
ties  are  computed  for  lattice  elements  with  the 
indices  n  and  m. 

According  to  Kerwin  (_1_1 ) ,  the  induced  velocity 
due  to  the  spanwise  vortices  and  sources  is 
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where  N  and  M  are  the  number  of  panels  over  the  chord 
and  span,  respectively. 

The  induced  velocity  due  to  the  chordwise 
vortices  and  the  separated  sheet  at  the  tip  is 
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where  Che  symbol  £  In  Che  second  summation  represents 
the  individual  chordwise  vortex  elemencs  due  to  the 
outer  end  of  the  n'th  spanwise  vortex. 

The  Induced  velocity  due  to  the  transition  wake 
is 
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where  N^m)  stands  for  Che  number  of  points  describing 
the  path  of  the  trailing  vortex  shed  from  the  inner 
end  of  rNmCsl  the  inner  radius  of  the  transition 
wake  should  be  exactly  the  same  as  the  hub  geometry. 
The  collection  point  of  the  separated  cip  vortex  con¬ 
tributes  the  (M+l)sc  transition  wake  element. 

The  induced  velocity  due  to  the  ultimate  tip 
vortex  is 


Nu-1 

Z 

n-1 


(5) 


where  Nu  represents  Che  number  of  points  which  form 
Che  piecewise  linear  approximation  to  Che  helical 
ultimate  wake.  The  first  poinc  of  the  helical  ulti- 
nace  wake  agrees  with  the  last  point  in  Che  tip 
element  of  the  transition  wake. 

Finally,  according  to  Wang  (10),  the  induced 
velocity  due  Co  the  hub  vortex  is 

_  NT  NT 

Vh  -  I  I  ri  ^ii  •  <*) 

1-1  J-l  J  J 

where  •  is  the  velocity  Induced  at  the  1th  control 
poinc  by  the  jth  panel.  NT  -  NR  *  NH  is  the  total  num¬ 
ber  of  panels.  Each  panel  consists  of  two  helical 
vortex  elemencs  and  two  vortex  ring  elemencs.  NH  is 
Che  number  of  helical  vortices  on  the  hub  between  two 
blades.  NR  -  NN1  +  NN  +  NN2  is  the  number  of  vortex 
rings,  where  NN1  is  the  number  of  vortex  rings  between 
the  blade  leading  edge  and  the  blade  trailing  edge,  and 
NN2  is  the  number  of  vortex  rings  between  the  blade 
trailing  edge  and  the  hub  apex.  The  total  Induced 
velocity  due  to  blade  and  wake  and  hub  vortices  is 

V“s  ♦  Te  +  Ttw  +*Vt  +  7h  .  (7) 

A  flow  chart  for  the  new  CR  lifting-surface  pro¬ 
gram  is  shown  in  Fig.  5.  The  lnpuc  Includes  the  blade 
characteristics  (such  ars  blade  and  hub  geometry,  and 
circulation  from  the  results  of  lifting-line  calcula¬ 
tions),  inflow  information,  wake  parameters,  and 
locations  at  which  field  point  velocities  Induced  by 
one  propeller  are  obtained,  Che  field  point  velocities 
on  the  ocher  propeller  are  automatically  calculated. 
These  field  point  velocities,  which  result  from  two 
sources  (one  from  blade  and  wake  vortices  and  the  other 
from  the  hub  vortex),  will  be  treated  as  additional 
onset  flows  when  Che  ocher  propeller  is  achieved,  Che 
program  will  automatically  compute  the  field  point 
velocities  Induced  on  the  first  propeller.  The  first 
propeller  will  be  redesigned  using  the  new  onset  flow. 
This  iterative  procedure  is  followed  until  Che  field 
point  velocities  induced  by  one  propeller  on  the  other 
show  no  significant  change. 

3.  Slipstream  Contraction 

As  flow  passes  through  a  propeller,  the  stream¬ 
lines  must  contract  because  the  propeller  Imparts  an 
axial  acceleration  to  the  flow  passing  through  the 
propeller.  The  contraction  of  the  slipstream  is  a 


well  known  phenomenon  from  the  evidence  of  theore¬ 
tical  analysis  as  well  as  experimental  measurements. 

It  is  very  difficult  to  accurately  predict  Che  tra¬ 
jectory  of  single  streamlines  analytically  due  to  the 
effects  of  nonlinear  and  three-dimensional  character¬ 
istics  of  the  flow.  This  is  especially  because  the 
vortex  sheets  tend  to  roll  up  in  the  slipstream.  Up 
to  this  point,  the  streamline  contraction  of  a  CR  set 
has  seldomly  been  Included  in  the  design  or  analysis 
of  CR  propellers.  Nagle  and  McMahon  (14)  measured  the 
flow  velocity  in  the  vicinity  of  a  sec  of  CR  research 
propellers  (propellers  4866  and  4867)  using  Laser 
Doppler  Velocimetry  (LDV).  The  propeller  sec,  designed 
to  operate  in  uniform  flow,  was  tested  in  the  DTRC 
24-in.  water  tunnel.  Axial,  radial,  and  tangential 
Induced  velocity  measurements  were  obtained  forward  of 
the  two  propellers,  between  the  propeller  set,  and 
downstream  of  the  aft  propeller,  with  the  propellers 
operating  at  the  design  condition. 


FIC.  5  FLOW  CHART  OF  NEW  CR  LIFTING- 
SURFACE  PROGRAM 

The  position  of  the  trailing  vortex  wake  geomecry 
is  traced  by  searching  for  the  aero  crossing  point  of 
the  tangential  velocity  (i.e.,  Vt-0)  versus  radius. 
Despite  the  large  number  of  data  points  which  were 
collected,  there  are  not  enough  data  points  to  deter¬ 
mine  the  complete  slipstream.  Daca  were  obtained  at 
two  different  locations  downstream  of  the  aft  propel¬ 
ler.  These  Include  experiments  17  through  21.  From 
the  data  collected  in  these  experiments,  one  slipstream 
behind  each  propeller  was  derived.  It  has  been  plotted 
in  Fig.  6.  The  contraction  angles  of  the  forward  and 
aft  propellers  are  15’  and  34°,  respectively.  It  is 
believed  chat  the  contraction  angle  of  the  aft 
propeller  is  larger  than  chat  of  the  forward  propeller 
because  the  Induced  velocity  effect  of  the  forward 
propeller  on  Che  aft  propeller  is  more  significant  than 
the  effect  of  the  aft  propeller  on  Che  forward  pro¬ 
peller. 

III.  SAMPLE  CALCULATION  FOR  A  CR  PROPELLER  DESIGN 

A  CR  propeller  set,  with  a  7  blade  forward  pro¬ 
peller  (D  -  16.5  feet)  and  5  blade  aft  propeller  (D  - 
15.826  feet),  has  been  designed  for  a  surface  ship. 

The  propellers  were  designed  for  a  speed  of  20  knots 
with  the  rotational  speeds  for  both  propellers  at  65 
rpm.  The  traditional  design  procedure  was  used  for 
this  design.  In  this  section,  a  comparison  between 
the  traditional  design  procedure  (TOP)  and  the  new 
design  procedure  (NDP)  t9  given. 
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FIG.  6  SLIPSTREAM  CC  TRACTION  OF  A  CR  PROPELLER  USING  LDV  MEASUREMENTS 
(This  figure  is  based  on  Nagle  and  McMahon  (14) ) 


For  Che  forward  propeller  design.  Fig.  7  shows 
the  radial  distribution  of  the  axial  velocity, 
uaA>p,  induced  by  the  aft  propeller  at  the  forward 
propeller  reference  line.  The  magnitude  of  UaA.F 
from  TOP  (lifting-line  model)  is  smaller  than  that 
from  NDP  (che  lifting-surface  model)  because  NOP 
Includes  the  thickness,  rake  and  skew  effects.  This 
indicates  that  the  torque  of  NDP  will  be  larger  than 
that  of  TOP.  However,  the  magnitude  of  uaAipis  very 
small  because  che  wake  induced  by  the  aft  propeller 
does  not  significantly  affect  the  forward  propeller. 


The  variation  of  the  radial  velocity,  UrA  p,  induced 
by  the  aft  propeller  at  the  forward  propefler  reference 
line  is  shown  in  Fig.  8.  The  magnitude  of  the  radial 
velocity  from  the  TDP  is  equal  to  zero,  while  that  of 
NDP  is  small  and  oriented  outward.  This  indicates 
that  the  amount  of  wake  contraction  for  the  forward 
propeller  will  be  reduced  due  to  the  effect  of  che 
aft  propeller.  Fig.  9  shows  that  the  tangential 
velocity,  utA  p, induced  by  the  aft  propeller  at  the 
forward  propeller  reference  line  Ls  zero  because  the 
wake  of  the  aft  propeller  can  not  affect  the  forward 
propeller. 


NON-DIMENSIONAL  RADIUS.  r/R 


FIG.  7  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  AXIAL  VELOCITY  INDUCED 
BY  AFT  PROPELi-ER  ON  FORWARD  PROPELLER 
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VELOCITY,  u, 


NON-OIMENSIONAL  RADIUS.  r/R 


FIG.  8  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  RADIAL  VELOCITY 
INDUCED  BY  AFT  PROPELLER  ON  FORWARD  PROPELLER 
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FIG.  9  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  TANGENTIAL  VELOCITY 
INDUCED  BY  AFT  PROPELLER  ON  FORWARD  PROPELLER 
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Since  the  differences  In  the  Induced  velocities 
predicted  by  TDP  and  NDP  are  very  small  for  the  for¬ 
ward  propeller,  the  pitch  and  the  camber  distributions 
should  be  similar.  The  radial  distribution  of  the 
pitch-diameter  ratio,  P/D,  is  affected  by  both  the 
Inflow,  the  propeller  induced  velocities,  and  by  the 
rotational  sneed.  3ecause  of  the  identical  rotational 
speeds  and  t.,.  low,  the  induced  velocities  dominate  the 
differences  in  the  pitch-diameter  ratio,  shown  in 
Fig.  10,  between  the  propellers  designed  using  TDP 
and  NDP.  The  magnitude  of  P/D  from  NDP  is  larger 
than  that  from  TDP  because  pfrom  NDP  la  larger. 

To  maintain  the  same  thrust,  the  camber  ratio,  f^/c, 
of  NDP  should  be  smaller  than  that  of  TDP  except  near 
the  hub.  This  is  because  NDP  Includes  hub  image 
effects  (see  Fig.  11).  The  thrust  coefficient, 
the  torque  coefficient,  ,  and  the  open  water  effi¬ 
ciency,  *0,  are  given  as  a  function  of  drag 
coefficient,!^,  in  Fig.  12.  nQ  from  NDP  is  smaller 
than  that  of  TDP  because  the  torque  of  NDP  is  larger 
than  that  of  TDP. 

As  for  the  aft  propeller  design.  Fig.  13  shows 
the  axial  velocity,  ua p  A, induced  by  the  forward  pro¬ 
peller  at  the  aft  propeller  reference  line.  The 
magnitude  of  ua  p  Afrom  TDP  is  smaller  than  that  from 
NDP  near  both  the  hub  and  the  tip  because  NDP  includes 
the  thickness,  rake  and  ske*.-.*  effects.  Thus,  it  i9 
expected  that  NDP  should  have  higher  torque  than  TDP. 
The  magnitude  of  ua p  A  is  higher  than  uaA  p because 
the  wake  of  the  forward  propeller  affects  the  aft 


propeller  more  than  that  of  the  aft  propeller  affects 
the  forward  propeller.  As  Fig.  14  shows,  the  radial 
velocity,  11^. p  ,  #  induced  by  the  forward  propeller  at 
the  aft  propeller  reference  line  has  zero  value  for 
TDP  because  the  lifting-line  model  does  not  account 
for  the  radial  Induced  velocity.  The  NDP  has  an 
inward  value  which  will  Increase  the  wake  contraction 
of  the  aft  propeller.  The  variation  of  the  tangencia 
velocity,  U£ p  A,  induced  by  the  forward  propeller  at 
the  aft  propeller  reference  line  is  shown  in  Fig.  15. 
The  magnitude  of  the  tangential  velocity  from  NDP  is 
higher  near  the  tip  and  lower  near  the  hub. 

The  pitch-diameter  ratios  from  the  designs 
developed  with  TDP  and  NDP  demonstrate  differences 
near  the  hub  but  not  the  tip,  as  shown  in  Fig.  lb. 

It  is  apparent  that  the  higher  pitch  from  NDP  near 
the  hub  is  caused  by  the  higher  uap  Aand  the  lower 
utF  A^rom  ^DP  c^an  fro®  TDP.  On  th’e  contrary,  the 
smafl  differences  in  the  pitch  near  the  tip  are  due 
to  Increases  in  both  i^p  A  and  p  A  when  computed 
by  NDP.  Fig.  17  shows  trte  radial  distribution  of 
the  camber,  fM/c.  It  can  be  seen  chat  the  magnitude 
of  camber  from  NDP  is  lower  than  that  from  TDP, 
expecially  near  the  tip.  The  reverse  camber  near 
the  tip  may  be  caused  by  the  sharp  decrease  of  the 
chord  distribution  near  that  region.  Fig.  18  shows 
ICj.,  Kq,  and  n0  versus  C^.  The  value  of  from  NDP 
is  smaller  than  that  from  TDP  because  the  thrust 
from  NDP  increases  but  the  torque  Increases  even 
more. 


FIG.  10  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  PITCH/DIAMETER 
RATIO  FOR  FORWARD  PROPELLER 
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'PEN  WATER  EFFICIENCY, , 
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RADIAL  INDUCED  VELOCITY,  u,  ,  A  AXIAL  INDUCED  VELOCITY,  u.  ,A 


FIG.  13  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  AXIAL  VELOCITY  INDUCED 
BY  FORWARD  PROPELLER  ON  AFT  PROPELLER 
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FIG.  14  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  RADIAL  VELOCITY  INDUCED 
BY  FORWARD  PROPELLER  ON  AFT  PROPELLER 
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PIT CH/CHAMETER  RATIO.  P/O  ■«  TANGENTIAL  INDUCED  VELOCITY,  u, 


.  15  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  TANGENTIAL  VELOCITY  INDUCED 
BY  FORWARD  PROPELLER  ON  AFT  PROPELLER 


FIG.  16  COMPARISON  OF  RADIAL  DISTRIBUTION  OF  PITCH/DIAMETER  RATIO 
FOR  AFT  PROPELLER 
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IV.  SUMMARY 

A  new  lifting-surface  computer  program  for  the 
detailed  deaigh  of  the  individual  propellers  of  a  CR 
set  has  been  developed  (see  Chen  and  Reed  ( 15 ) ) .  In 
the  traditional  design  procedure,  Che  lifting-surface 
program  computed  the  induced  velocities  from  one 
propeller  on  the  other  from  lifting-line  calculations. 
The  new  design  program  computes  the  induced  velocities 
from  lifting-surface  cheory.  Additionally,  the  hub 
portion  of  the  program  has  been  modified  to  take  into 
account  the  velocities  Induced  by  one  propeller  on  the 
hub  of  the  other  propeller.  The  shape  and  distribution 
of  the  wakes  shed  from  the  two  propellers  are  adjusted 
based  on  LDV  measurements  of  Induced  velocities. 
Nevertheless,  the  adjustment  is  very  crude  due  to  a 
limited  amount  of  data. 

A  sample  calculation  which  shows  the  comparison 
between  the  traditional  and  the  new  design  procedures 
was  shown.  The  results  from  the  two  design  procedures 
look  fairly  consistent  with  each  other  except  near  the 
hub  region.  This  consistency  is  because  both 
procedures  input  che  same  circulation  distribution. 
However,  adding  hub  effects  at  the  lifting-line  step 
can  produce  a  much  different  circulation  distribution 
and  the  results  should  be  quite  different.  The  present 
study  indicates  that  the  use  of  Interaction  velocities 
calculated  from  che  lifting-line  model  is  sufficiently 
accurate  for  design  use  when  the  propeller  loadings  are 
moderate,  and  che  propellers  are  not  too  close  to  each 
other. 

V.  CONCLUSION'S  AND  RECOMMENDATIONS 

The  conclusions  from  this  study  are  as  follows. 

1.  The  new  lifting-surface  design  procedure,  which 
incl"d»s  automatic  Induced  velocity  computation,  hub 
effects,  and  empirical  adjustment  for  wake  contraction, 
is  a  better  propeller  design  method  than  previously 
used . 

2.  The  capability  of  thia  new  lifting-surface 
program  Is  not  restricted  co  CR  propellers  but  can  be 
used  to  design  ocher  compound  marine  propulsors  such 
as  vane  wheel  propulsors. 

There  are  several  recommendations  for  further 
study  of  this  topic. 

1.  The  effects  of  propeller  load  distribution  and 
spacing  need  to  be  further  investigated  to  determine 
when  the  use  of  lifting-surface  theory  is  mandated  for 
predicting  che  mutual  interaction  velocities. 

2.  Detailed  LDV  measurements  of  the  flow  velocity 
in  the  vicinity  of  CR  propellers  need  to  be  performed 
to  get  further  information  on  the  shape  and  distribu¬ 
tion  of  the  wakes  shed  from  the  two  propellers. 

3.  An  unsteady  force  calculation  design  method 
for  CR  propellers  needs  to  be  developed. 
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2  HQS  COGARD 

1  US  COAST  GUARD  (G-ENE-4A) 

1  LC/SCI  &  Tech  Div 

2  NASA  STIF 

1  DIR  RES 

1  NSF  ENGR  DIV/Lib 

1  DOT  Lib 

2  U  CAL  BERKELEY/DEPT  NAME 

1  NAME  Lib 
1  Webster 

1  U  MISSISSIPPI  DEPT  OF  M.E. 

1  Fox 


NAVSEA  3  CIT 


5 

SEA 

05R 

1 

AERO  Lib 

1 

SEA 

55 

1 

Acosta 

3 

SEA 

55N 

1 

Wu 

1 

SEA 

55W 

3 

SEA 

56D 

4  U  IOWA 

1 

SEA 

56X 

1 

Lib 

3 

SEA 

56X1 

1 

I HR/ Kennedy 

1 

SEA 

56X2 

1 

IHR/Stern 

3 

SEA 

56X4 

1 

1HR/ Patel 

1 

SEA 

56X5 

1 

SEA 

56XP 

1 

PMS- 

■383 

1 

PMS- 

•392 
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INITIAL  DISTRIBUTION  (Continued) 


Copies 

Copies 

3 

U  MICHIGAN/DEPT  NAME 

1 

WHO I  OCEAN  ENGR  DEPT 

1  NAME  Lib 

1  Parsons 

1 

WPI  ALDF.N  HYDR  LAB  Lib 

1  Vorus 

1 

ASME/RES  COMM  INFO 

3 

MIT 

1  BARKER  ENGR  Lib 

1 

ASNE 

2  OCEAN  ENGR/Kerwin 

1 

SNAME/Tech  Lib 

3 

STATE  U  MARITIME  COLL 

1  ARL  Lib 

1 

AERO  JET-GENERAL/Beckwith 

1  ENGR  DEPT 

1  INST  MATH  SCI 

1 

ALLIS  CHALMERS,  YORK,  PA 

4 

PENN  STATE  U  ARL 

1 

AVCO  LYCOMING 

1  Lib 

1  Henderson 

1  Gearhart 

1 

BAKER  MANUFACTURING 

1  Thompson 

1 

BIRD-JOHNSON  CO/Norton 

1 

BOEING  ADV  AMR  SYS  D1V 

1 

DOUGLAS  AIRCRAFT/Lib 

1 

BB&N/Jackson 

2 

EXXON  RES  DIV 

1  Lib 

1 

BREWER  ENGR  LAB 

1  Fitzgerald 

1 

CAMBRIDGE  ACOUS/Junger 

1 

FRIEDE  &  GOLDMAN/Miche 1 

1 

CALSPAN,  INC/Ritter 

1 

GENERAL  DYNAMICS,  EB/Boatwr ight 

1 

STANFORD  U/Ashley 

1 

GIBBS  A  COX/Lib 

1 

STANFORD  RES  INST  Lib 

1 

ROSENBLATT  &  SON/Lib 

2 

SIT  DAVIDSON  LAB 

1 

GRUMMAN  AEROSPACE/Car 1 

1  Lib 

I  Savitski 

1 

TRACOR  HYDRONAUTICS/Lib 

1 

TEXAS  U  ARL  Lib 

1 

INGALLS  SHIPBUILDING 

1 

UTAH  STATE  U/Jeppson 

1 

INST  FOR  DEFENSE  ANAL 

2 

VPI/DEPT  AERO  &  OCEAN  ENGR 

1 

ITEK  VIDYA 

1  Schetz 

1  Kaplan 

1 

LIPS  Duran/Kress 

2 

WEBB  INST 

1 

LITTLETON  R  &  ENGR  CORP/Reed 

1  Ward 

1  Hadler 

1 

LITTON  INDUSTRIES 
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INITIAL  DISTRIBUTION  (Continued) 


Copies 

Copies 

Code 

Name 

1 

152 

Lin 

1 

LOCKHEED,  SUNNYVALE/Waid 

1 

1521 

Day 

2 

MCDONNEL  DOUGLAS,  LONG  BEACH 

1 

1521 

Karaf iath 

1  Cebeci 

1 

1521 

Hurwitz 

1  Hess 

1 

1522 

Remmers 

1 

MARITECH,  INC/Vassilopoulis 

1 

1522 

Wilson 

2 

HRA,  INC 

1 

154 

McCarthy 

1  Cox 

1  Scherer 

1 

154.1 

Y  im 

1 

NIELSEN  ENGR/Spangler 

1 

1542 

Huang 

1 

NKF  ASSOCIATES/ Noonan 

1 

1544 

Reed 

20 

1544 

Chen 

1 

NAR  SPACE/Ujihara 

1 

156 

Cies lowski 

2 

ATLANTIC  APPLIED  RESEARCH 

1  B  rown 

1 

172 

Krenzke 

1  Greeley 

1 

1720.6 

Rockwell 

1 

PROPULSION  DYNAMICS,  INC 

1 

19 

Sevik 

1 

SPERRY  SYS  MGMT  Lib/Shapiro 

1 

1901 

Strasberg 

1 

TETRA  TECH  PASADENA/Furuya 

1 

1905 

Blake 

1 

UA  HAMILTON  STANDARD /C orne 1 1 

1 

1906 

Biancardi 

CENTER  DISTRIBUTION 

1 

1942 

Mathews 

Copies 

Code  Name 

1 

1962 

Kilcullen 

1 

0120 

10 

522.4 

Reports  Control 

1 

12 

1 

522.1 

TIC  (C) 

1 

112.1  Nakonechny 

1 

522.2 

TIC  (A) 

1 

15  W.R.  Morgan 

1 

1506  Walden 

1 

1508  Boswell 
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